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We suggest to view ultracold atoms in a time-periodically shifted optical lattice as a "dressed mat- 
ter wave", analogous to a dressed atom in an electromagnetic field. A possible effect lending support 
to this concept is a transition of ultracold bosonic atoms from a superfluid to a Mott-insulating state 
in response to appropriate "dressing" achieved through time-periodic lattice modulation. In order to 
observe this effect in a laboratory experiment, one has to identify conditions allowing for effectively 
adiabatic motion of a many-body Floquet state. 



I. INTRODUCTION 

The "dressed- atom picture" provides a transparent ap- 
proach to the dynamics of atoms and molecules in strong 
electromagnetic fields [1] . In essence, the "dressing" pro- 
vided by the field may equip the atom or molecule with 
properties quite different from those of a "bare" one. 
A hallmark example along these lines is given by the 
modification of Zeeman hyperfine spectra of atoms in- 
teracting with strong radiofrequency fields: As reviewed 
in Sec. IIIII below, in the presence of the dressing field 
the bare atomic ^-factors become multiplied by a certain 
Bessel function, the argument of which is proportional to 
the strength of the field, and inversely proportional to its 
frequency. Thus, experimentally recorded hyperfine spec- 
tra depend sensitively on these parameters y. Closely 
related phenomena have been observed, for instance, in 
radiatively assisted collisions of Rydberg atoms Q . 

The appearance of a Bessel function in response to 
time-periodic forcing is typical for quantum systems that 
can be viewed as a set of interacting nearest neighbours, 
such as the angular momentum substates in the case 
of the ^-factor modification. Another striking example 
has emerged only recently: In experiments with ultra- 
cold atoms in time-periodically shifted optical lattices, 
the hopping matrix elements which quantify the mag- 
nitude of the tunnelling contact between states located 
at adjacent lattice sites differ from those of a bare lat- 
tice system again by a Bessel function, and thus can be 
tuned by adjusting the amplitude or the frequency of the 
lattice modulation 0, . This finding now suggests an 
interesting question: The (7-factor experiment [3] with 
single atoms has been instrumental for establishingthe 
dressed-atom picture; could the recent experiments 0, Q 
with Bose-Einstein condensates lead to a similar picture 
of "dressed matter waves"? In other words, can one ex- 
ploit time-periodic forcing for endowing a macroscopic 
matter wave with properties it does not have when the 
forcing is absent? 

In order to expand on this question, we proceed as 
follows: We first recall in Sec. [Til the physics underlying 



the Bessel-function modification occurring in the dressed- 
atom picture, using the example of a two-level system 
interacting with a quantised radiation mode. However, 
when dealing with cold atoms in time-periodically shifted 
optical lattices, it is certainly reasonable to describe the 
time-periodic lattice modulation in terms of an exter- 
nal classical force. Hence, we discuss in Sec. IIIII how 
the Bessel function appears in that context, employing 
the Floquet formalism. For fully exploiting the possibil- 
ities of control opened up by time-periodic forcing, adi- 
abatic response to slowly changing parameters plays an 
important role. Therefore, we briefly point out in Sec. lIVI 
how the adiabatic principle works for Floquet states. In 
Sec.|V]we put all these pieces together and argue that a 
system of interacting ultracold bosonic atoms in a time- 
periodically modulated optical lattice can change its state 
from superfiuid to Mott-insulator-like and back, if the 
modulation is switched on and off in an adiabatic man- 
ner, and the parameters are chosen judiciously @, 0iS]- 
This scenario, still to be confirmed experimentally, could 
lend significant support to the notion of dressed matter 
waves: The dressing achieved through the lattice modu- 
lation determines the state of the system. We finally sum 
up our conclusions in Sec. IVII 



II. THE DRESSED TWO-LEVEL SYSTEM 

We start by studying a "two-level atom" interacting 
with both a static external field and a single mode of a 
quantised radiation field. Its dynamics are governed by 
the Hamiltonian 



where 
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(1) 



(2) 



specifies the unperturbed system with energy eigenvalues 
±J/2, 
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models a radiation mode with frequency w in terms of a 
harmonic oscillator with dimensionless oscillator coordi- 
nate z, and 



Hint = i (i^O + IZ) ' ^ 





-1 



(4) 



describes the static field of strength Kq and the coupling 
to the radiation mode, with a strength specified by a con- 
stant 7. When J = 0, so that the two "atomic" levels 
are degenerate, this Hamiltonian ([T|) obviously is diago- 
nalised by the shifted harmonic-oscillator states 



with energies 



ipn{z + -1 1 21110) \ 

) 


ipn{z - ^/2huj) 



(5) 



(6) 



the functions (fniz) denote the familiar eigenfunctions of 
the harmonic oscillator, 

^„(z) = (V¥2" ^^^(^) exp(-zV2) • (7) 

The static field splits the two atomic levels by the amount 
AE — Kq. When this splitting is matched by an integer 
number of photons, that is, when 



Ko = {i-n)huj , 



(8) 



the unperturbed states are pairwise degenerate, — 

E^^l = 0. We now study the removal of this degeneracy 
between an "n-photon state" and an "^-photon state" for 
nonzero J, assuming J <C 

Degenerate-state perturbation theory requires to eval- 
uate the matrix elements of the "perturbation" iJat in the 
basis ([5]), and, hence, to compute the overlap integrals 



M= dzipn{z + j/2huj)ipi{z-^/2huj) . (9) 

^ — 00 

With the help of the expansion 

fc=0 ^ 

this integral can be calculated exactly, yielding 



(10) 



M = exp(-a/2) a 



k=0 



where we have introduced the dimensionless parameter 

'^='{2ty- 



kl (n~k)l {£-k)\ 

(11) 
ter 

(12) 



Assuming £ > n, and employing the generalised Laguerre 
polynomials [9|] 



fc=0 



kl {n- ky.{e^n + ky. 



(13) 



this expression pT|) takes the form 



M = (-l)^-"exp(-a/2)a(^-")/^ V ^"""^^"^ ' 

(14) 

Now we are interested in the limiting case of almost clas- 
sical fields containing a very large number of photons. In 
order to maintain the resonance condition ([8]), we keep 
the integer v = I — n fixed while letting n and £ tend to 
infinity. In that limit, one has Q 

L(f-")(x) "-°°^"='' lle-/2 (^^±^^Y'~^^'^ 
X j£_„(v/2(n + ^+l)2;) , (15) 

where i,y{x) denotes a Bessel function of integer order v. 
Hence, one finally obtains 



M 



(-l)^-"J,_„(V2(n + ^-f l)a) (16) 



in that same limit. 

It remains to interpret the peculiar-looking argument 
of the Bessel function. When placing the field oscillator 
into a coherent state, its amplitude zq is determined by 
energy considerations: Since, as expressed by the reso- 
nance condition ([5]) , the field is exchanging t — n photons 
with the atom, the average field energy is the arithmetic 
mean of the energy of an n-photon state and that of an 
^-photon state. This gives 



1 



which, in view of the definition (jl2p . implies 



7^0 
Tiu) 



Da 



(17) 



(18) 



The energy eigenvalues which have been degenerate for 



J 0, i.e., E^°]^ = El'!_, now are shifted by ±MJ/2 for 
nonzero J. For energies sufficiently high to validate the 
preceding reasoning, the spectrum of the Hamiltonian ([1]) 
thus consists of a sequence of doublets spht by \MJ\, 
the doublet centers being separated by fiLO. Putting all 
things together, this means that the "atom" iJat, when 
"dressed" by the interaction H^^i + -^rad under the con- 
ditions detailed above, behaves like a noninteracting sys- 
tem ^ with a modified level splitting determined by the 
effective J-parameter 



AO) 



j.» = (-irj.(S)j. 



(19) 
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This is the lesson to be learned from the present two- 
level example: When the driving field can be considered 
classical, resonant forcing effectuates a modification of 
the "atomic" level splitting such that the unperturbed 
splitting is multiplied by a Bcssel function with an argu- 
ment proportional to the driving amplitude divided by 
huj. The order v of this Bessel function is determined 
by the number of photons in resonance with the transi- 
tion, according to Eq. ([8]). In particular, when there is 
no static field, one has 1^ = and thus recovers the mod- 
ification of the splitting by a Jq Bessel function which 
also underlies, for instance, the coherent destruction of 
tunnelling of a sin gle p article in driven symmetric double 
well potential [lol llll|. Quite recently, this phenomenon 
has been observed with cold atoms in periodic double- 
well potentials [T2|] . 



III. ELEMENTS OF FLOQUET THEORY 

In order to avoid the consideration of a quantised field 
and to start with a classical driving force right away, we 
now treat the explicitly time-dependent Hamiltonian 



H{t)=Ho + Hi{t) , 



(20) 



where the time-independent system Hq corresponds to a 
spin 1 in a magnetic field B oriented in the x-direction. 




(21) 



with /xb denoting the Bohr magneton, and gi the Lande 
g-factor. The external forcing is given as an additional 
static magnetic field Bq and an oscillating field with am- 
plitude Bi^ and frequency uj, both directed along the z- 
axis, 

/ 1 

ffi(<) =5iMB[Bo + B^cos(wt)] I . (22) 

\0 0-1 

Apart from the fact that here the forcing is truly classi- 
cal, this system closely resembles the previous two-level 
example ((T]): The unperturbed system is characterised 
by "nearest-neighbour coupling" , while the forcing is di- 
agonal. Now the Hamiltonian ([20]) depends periodically 
on time, 

H{t) = H{t + T), (23) 
with period T = 2Tr/(jj. Hence, the Schrodinger equation 
d 



ih-m)) ^ H{t)\m) 



(24) 



has Floquet-type solutions [ll [13, d, [H 

l^a(O) = W{t)) exp(-ie„V?i) , (25) 



where the functions \ua{t)) inherit the periodic nature of 
H{t), 



\uo.(t)) - \u^{t + T)) 



(26) 



These functions, together with the corresponding 
quasienergies Eq,, are obtained as solutions to the eigen- 
value problem 

{H{t) - ihdt) \u^{t))) = e„ |m,(0)) , (27) 

which is defined in an extended Hilbert space TL ®T oi 
T-periodic functions [171] in which the time t is regarded 
as a coordinate and which, therefore, is equipped with 
the scalar product 



/ di(.|.). 



(28) 



combining the standard scalar product ( • | • ) for the sys- 
tem's original Hilbert space TL with time-averaging. We 
stick to the convention of writing \ua{t)) for a Floquet 
function viewed in 7i, but \ua{t))) when that same func- 
tion is regarded as an element of the extended space 

n®T. 

There is one issue implied by the Floquet formalism 
which requires particular attention. Namely, if |u(-„ Q)(t))) 
solves 

{H{t) - indt)\u(^^o){t))) = £(„^o)|w(n,o)(i))) (29) 
with quasienergy £(„.o), then 

\u{n,m){t))) = |w(n,0)(i))) exp(iTOWt) (30) 

solves 

[H{t) - \Tldt)\u(n,r,i){t))) = £(n,m)\u(n,rn){t))) (31) 

with quasienergy 

e(n,m) = £{nfi) + mhuj , (32) 

where m is any (positive or negative) integer. Hence, the 
quasienergy spectrum repeats itself periodically on the 
energy-axis; each "Brillouin zone" of width huj contains 
one respresentative, labelled by m, of the class of eigen- 
values belonging to the Floquet state labelled by n. But 
when following the evolution of a wave function \'ip{t)) 
in the physical Hilbert space TL, only one representative 
from each class is needed, giving an expansion of the form 

IV'(0> = X!*^" l"(n,o)(i)) exp(-ie(„^o)V^) (33) 

n 

with time-independent coefficients c„. 

In order to apply this lore to the spin-l-system (|20p . 
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we observe that the Floquet basis states 



|w(0.m)(i)» 
\u{-,m)it))) 



I exp ( — — — sin(u;i) + imujt 





1 I exp(imti)t) 




I exp ( — — sin(wi) + imujt 



hui 



(34) 



diagonahse the quasienergy operator Hi{t)—\hdt which is 
obtained when there is no field B, so that the eigenstates 
of the three-level Hamiltonian ([2T|) are degenerate; the 
"unperturbed" quasienergies express the Zeeman split- 
ting caused by the other static field Bq: 



(n = 0,±11 



If the oscillating field is resonant in the sense that 
then the Floquet functions 

= l"(+,o)W)) 

\u2[t))) = |?/(0,.)(t)» 
\Mt))) = l"(-,2.)W» 



(35) 



(36) 



(37) 



correspond to the same quasienergy, and thus are de- 
generate. The removal of this degeneracy for nonvan- 
ishing B once again is assessed by degenerate-state per- 
turbation theory, assuming gi/ie-B ^ fiw. In contrast 
to Sec. nil now the calculation proceeds in the extended 
Hilbert space HC^T, and thus invokes the computation 
of matrix elements {{uj\Ho\uk)) with respect to the scalar 
product (|28l) . But this is what makes the mathematics 
quite simple: Using the identity 



^ e"^-*Jfc(z) 



(38) 



00 



for expanding the unperturbed Floquet functions, time 
averaging according to the definition (pS)) serves to filter 
out one particular term from the sum, determined by the 
resonance condition ([36]). Thus, one immediately obtains 



{{u,\Ho\uk)) = {- 



B, 



V flUJ 



1 



V2 



J,fc±l 



giving the quasienergies 



f aim 



B^ 



\ flLO 



(39) 

mhuj . (40) 



Hence, the effect of the forcing (|22p on the system (pT|) 
is described by replacing the bare 5-factor gi by the ef- 
fective substitute 



ffeff = (-l)''Ji. 



/ .91 Mb -Be 



91 ■ 



(41) 



Evidently, the line of reasoning adopted in this section 
to treat classical forcing parallels the arguments given 
in Sec. inifor a system interacting with a quantised field. 
But here the argument is considerably more direct, avoid- 
ing the analysis referring to "large photon numbers" . The 
price to pay for this simplification is a quasienergy spec- 
trum which is strictly ?ia;-periodic und thus unbounded 
from below, whereas the exact quantum mechanical en- 
ergy spectrum becomes approximately ?ia;-periodic only 
for sufficiently high quantum numbers. Nonetheless, for 
systems subjected to time-periodic classical forcing the 
Floquet picture combines great conceptual clarity with a 
fairly succinct computational approach. 



IV. ADIABATIC FOLLOWING OF FLOQUET 
STATES 



One additional piece of input is required before we can 
treat ultracold atoms in a periodically shifted optical lat- 
tice, namely, the adiabatic response of Floquet states to 
slowly changing parameters. In order to make the point, 
let us briefly recapitulate the standard adiabatic theorem 
of quantum mechanics [ISj]: The task is solve a time- 
dependent Schrodinger equation 



(42) 



with a Hamiltonian H^^^^ depending on a parameter P{t) 
which changes slowly in time. The strategy then is to 
"freeze" that parameter in a first step, and to consider 
the family of eigenvalue problems 



(43) 



for each relevant, fixed value of P. Let us stipulate that 
the phases of the instantaneous eigenstates \ be cho- 
sen such that 







(44) 



If then the system initially, at time t = 0, is prepared in 
a particular eigenstate. 



m = 0)) = i<^^(*=°)) , 



(45) 



and P is allowed to vary sufficiently slowly, an approxi- 
mate solution to the Schrodinger equation is given 
by 



|^(0) = |^r^*))exp(-i^dt'<(*': 



(46) 



provided the parameter variation proceeds smoothly, and 
|iy9^) is separated for all P by an energy gap from the 



5 



other states. Hence, the system stays in the state contin- 
uously connected to the one it was originaUy prepared in, 
and acquires a "dynamical" phase determined by an in- 
tegral over the instantaneous energy eigenvalues encoun- 
tered during its evolution. We remark that it might not 
be possible to satisfy the phase-fixing condition (j44|) glob- 
ally if there is more than one time-dependent parameter; 
this fact then forces one to explicitly introduce Berry's 
geometrical phase [lOj . 

When trying to transfer this adiabatic theorem to 
systems with a Hamiltonian H^^*^{t) which would de- 
pend periodically on time if the parameter P where 
fixed, H^{t) = H^{t + T), but which actually exhibits 
an additional "slow" time-dependence of P, one faces 
a problem: If one simply "stopped the time" in order 
to define an instantaneous Hamiltonian, one would not 
only freeze the parameter P, but also loose the peri- 
odic time-dependence. However, it appears much more 
natural to freeze only P, and to maintain the periodic 
time-dependence on the level of the instantaneous eigen- 
value problems. The way to do so, as formulated in 
Refs. [13, HH, includes a detour to the extended Hilbert 
space 7i (g) T introduced in the previous section: Instead 
of starting from the actual Schrodinger equation 



(47) 



one first distinguishes two different time variables, a vari- 
able T for the slow, parametric time dependence and a 
variable t for the fast, oscillating one, and then considers 
the evolution equation 

ihdr\^iT,t))) - (i/^(-)(i) -iMt)|*(r,0)) (48) 

in §5 T. If this equation can be solved, one returns to 
the desired wave function \tp(t)) evolving in the system's 
true Hilbert space H by equating t and t: One has 



(49) 



since 



mt\m) - ihdr\-^{T,t))) 



{HP(^\t)~ihdt)\^{T,t))) 



(50) 



On the level of the extended evolution equation ([481) , one 
can now freeze P by stopping solely the time r, while 
leaving the other time t unaffected. This then defines 
the instantaneous eigenvalue problems in terms of the 
operators appearing on the right-hand side of Eq. (|48)) . 



(51) 



Since, by construction, this problem lives in 7i (g) T, it is 
exactly the quasienergy problem formulated in Eq. (|27p . 



The remaining reasoning follows the standard route: We 
fix the phases of the instantaneous eigenstates by requir- 
ing 

(«|5pO)=0, (52) 
and start at time r = with the initial condition 



|vI/(r = 0,t))) = K(^=°)(i))) 



(53) 



Then 



|vl/(r,t))) = |u^M(t)))exp(^-i dr'e^(^')j (54) 

is an adiabatic solution to the extended evolution equa- 
tion (UHl) , provided the propositions of the adiabatic theo- 
rem can be met, and returning to Ti. according to Eq, 
gives 



(i))exp(^-i^di'£r<*') 



(55) 



as an approximate solution to the original Schrodinger 
equation (|47p . In short, for adiabatic quantum trans- 
port in periodically time-dependent systems with slowly 
changing parameters the Floquet states adopt a role 
which is completely analogous to that played by energy 
eigenstates in conventional situations described by an 
equation of the type (I42p . The strategy of "lifting" the 
Schrodinger equation (|42ll to the extended space (g) T, 
applying standard techniques there, and then projecting 
back to Ti. is useful not only for understanding the struc- 
ture of the problem, but also for detailed computations 
of non-adiabatic corrections [2^ . 

There is, however, a big caveat. As remarked above, 
the standard adiabatic theorem demands that the adi- 
abatically transported state be separated by an energy 
gap from all other states. Accordingly, when trans- 
ferring this theorem to H (E) T, one requires that the 
adiabatically transported Floquet state be separated 
in quasienergy from the other ones. But since one 
quasienergy-representative from each state falls into each 
quasienergy Brillouin zone, this condition is almost im- 
possible to satisfy when there is a large number of states. 
One then expects a multitude of near-degeneracies "mod- 
ulo huj" , reflecting a dense set of multiphoton resonances. 
In such a situation, it appears unlikely that an adiabatic 
limit exists j23| . However, it appears equally plausibe 
that, if one does not consider the fictitious limit of a 
parameter variation proceeding "infinitely slowly" , but 
instead specifies that the variation takes place within 
a finite time interval, most of these resonances are not 
"seen" long enough by the system to become active. 
Then effectively adiabatic motion is possible, if major 
resonances can be avoided. Although it might be hard 
to formulate this somewhat vague notion in a mathemat- 
ically precise manner in the general case, the emerging 
adiabatic principle (not theorem) for Floquet states can 
provide intuitively clear guidelines for understanding the 
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evolution of periodically driven systems in well-designed 
particular cases. The following discussion of the driven 
Bose-Hubbard model exemplifies that the occurrence of 
effectively adiabatic motion, or its destruction by active 
resonances, depends on the choice of the frequency. 

V. THE DRIVEN BOSE-HUBBARD MODEL 

The Bose-Hubbard model, as devised by Fisher et 
al. [l^l, describes Bose particles on a lattice. There ex- 
ists a tunnelling contact between neigbouring sites, with 
a strength specified by a hopping matrix element J; each 
pair of particles occupying the same site increases the 
energy of the system by an amount U due to repulsion. 
Thus, for the case of a one-dimensional (Id) lattice with 
M sites the many-body Hamiltonian reads 

Ha = -J {^y^i+i + ^^l+iW) + ^Y^eif^i-^) . 
e=i 1=1 

where b\ (b^) is the creation (annihilation) operator for a 
Bose particle at the ^th lattice site, obeying [b^, b\] — S^^k, 
and "hi = b\b^ gives the number of particles on that site. 
Assuming that there are N particles in total, and that 
the filling factor n = N/M is integer, the system's ground 
state undergoes a significant change when the dimension- 
less control parameter U/J is varied: In the interaction- 
free limit U/J ^ it corresponds to a superfluid, given 
by a Bose-Einstein condensate with all particles occupy- 
ing the lowest Bloch state, 

1 / 1 "A" 

1°). (5^) 

where |0) is the vacuum state. In the opposite limit of 
vanishing tunnelling contact, U/J — > oo, the individual 
sites are isolated, so that the systems adopts the Mott- 
insulating ground state 

m^Yl^-jl^io). (58) 

1=1 

When the lattice is infinitely large, that is, for M — > oo 
and N ^ oo while keeping n — N/M constant at an 
integer value, a sharp transition between the superfluid 
and the Mott-insulating regime occurs at a critical value 
{U/ J)c, accompanied by the emergence of a finite en- 
ergy gap. For a Id lattice with filling factor n — 1, 
one finds {U/J)c ^ 3A The Bose-Hubbard model 

has received considerable attention recently, since it can 
be realised with ultracold atoms in d-dimensional optical 
lattices {d — 1,2,3) [1^, allowing one to investigate the 
superfluid-to-Mott insulator quantum phase transition in 
great detail in the laboratory [13, [H, 0, [30| . 

We view ultracold Bose particles in an optical lattice 
as prime candidates for exploring the concept of dressed 



matter waves. Namely, atoms in a Id lattice can be sub- 
jected to a time-periodic lattice modulation, to the effect 
that a term of the form 

M 

Hi {t) = [Kn + cos(wt)] ^ thi (59) 

£=1 

is added to the system ([56| . Here denotes the ampli- 
tude of a drive with angular frequency typically on the 
order of one to a few kilohertz while Kq corresponds 
to a static lattice tilt [sl]; the extension to lattices with 
higher dimension is straightforward. If there were no in- 
teraction between the particles, that is, for U/J = 0, the 
total Hamiltonian H{t) — Hq + Hi{t) would be identical 
in form to the systems studied in Sees. [Til and IIIII One 
faces nearest-neighbour coupling combined with homoge- 
neous site-diagonal forcing. Hence, when the resonance 
condition 

Ko = i^huj (60) 

corresponding to the previous equations ^ and ([55]) is 
satisfied, so that the energy of "photons" matches the 
energy shift induced by the static tilt between adjacent 
sites, one can adapt the results (fTO]) and (|4T|) : Under such 
conditions, the driven system behaves approximately like 
an undriven one with the modified hopping matrix ele- 
ment 

J,ff^i~ir3AKu/hco)J. (61) 

More careful analysis @, 0| shows that this expression 
remains valid even for nonzero U at least in the high- 
frequency regime where hco 3> U and huj 3> J, thus 
including the strong-coupling case U/J ^ 1. The ex- 
perimental data available so far give clear evidence of 
this modification (|6ip both for v — 0, when there is no 
static tilt [4], and for = 1,2, when one has "photon"- 
assisted tunnelling This finding now directly leads 
to a further consequence: The ratio U/J governs the 
superfluid-to-Mott insulator transition in the bare Bose- 
Hubbard model; this control parameter has to be re- 
placed by U/JcB in the presence of resonant forcing. 
Since Joff depends significantly on the amplitude K^^, it 
should be possible to cross the border between the su- 
perfluid and the insulator regime by varying that ampli- 
tude @, 0]- However, the notion of a "superfluid" or a 
"Mott insulator" refers to the ground state of the Bose- 
Hubbard model, so that it becomes necessary to guide 
the ground state of the undriven system Hq into the ef- 
fective ground state of the driven system Hq + Hi{t). 
That "effective ground state", of course, is the Floquet 
state which originates from the ground state of Hq when 
the drive is turned on, so that the adiabatic principle dis- 
cussed in Sec. IIVI comes into play: In an experiment aim- 
ing at a demonstration of a superfluid-to-Mott insulator 
transition induced by time-periodic forcing, the driving 
amplitude should be turned on smoothly, such that the 
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FIG. 1: Quasimomentum distribution computed according 
to Eq. (|62[) for a driven Bose-Hubbard model with N = 7 
particles on M = 7 sites. The interaction strength is U/ J — 3, 
the driving frequency is Tiu/J — 16. Between t — and 
ti = 100 T the amplitude K^j (t) is increased linearly from 
zero to -ft'niax = 2.4 hui, then kept constant until t2 = 200 T, 
and finally ramped linearly back to zero between t2 and ta = 
300 T. Here the system is able to follow adiabatically. 



system's wave function can follow the Floquet state con- 
nected to the unperturbed ground state. But since adia- 
batic following in a periodically forced many-level system 
with a "dense" quasienergy spectrum is endangered by a 
host of multiphoton-like resonances, the precise choice of 
the protocol is not trivial: The frequency has to be chosen 
such that major resonances are avoided, while the ampli- 
tude has to vary sufficiently slowly in order to allow for 
a reasonable degree of adiabaticity, but still sufficiently 
fast in order to pass minor resonances before they become 
active. 

We illustrate these deliberations by numerical calcula- 
tions for small systems with N — 7 particles on M — 7 
lattice sites. The initial A^-body wave function \ip{0)) 
at time t = is chosen as the ground state of the Bose- 
Hubbard model ((56|) with U/J — 3, falling into the super- 
fluid regime when the system is sufficiently large. Here we 
restrict ourselves to i^o = 0, i.e., to = 0; related stud- 
ies for 1/ = 1,2 are documented in Rcf. [7|. After select- 
ing some frequency uj, and thus specifying the time scale 
T = 2tt/lo, the time-dependent force is turned on with an 
amplitude K^{t) which rises linearly between t = and 

= 100 T from K^/huj = to K^^Jhui = 2.4. The lat- 
ter value lies close to the first zero of Jq, and thus gives 
a quite large ratio U/ Jcs, which should place the sys- 
tem far into the Mott-like regime. Then the amplitude 
is kept constant at iiTmax between ti and ^2 = 200 T, 
and finally ramped linearly back to zero between t2 and 
ts = 300 T. The corresponding A^-body wave function 
\ip(t)) is computed by plain direct solution of the time- 
dependent Schrodinger equation, not taking any recourse 



FIG. 2: As Fig. [T] but for a lower frequency Tiuj/J = 12. 
In this case, adiabatic following is destroyed by Landau- 
Zener transitions at the avoided quasienergy crossings visible 
in Fig. H 



at all to Floquet theory. From that wave function, the 
single-particle quasimomentum distribution p{p, t) is ob- 
tained according to 



h/a 



(62) 



and recorded at integer multiples of T [1, 0] . This mo- 
mentum distribution is sharply peaked, due to (quasi) 
long-range phase coherence, in the superfluid phase, but 
apparently structureless in the Mott regime. 

Figure [1] shows the results ior hui / J = 16: Initially one 
finds a strongly peaked distribution, as expected for a 
superfluid-like state, which then becomes practically flat 
at t = ti. However, after the driving amplitude K^^it) 
has been switched off again at t = t^, the sharp pattern 
reappears: This fact clearly signals that the flat distri- 
bution between ti and t2 is not due to loss of coherence 
resulting from uncontrolled excitations, but rather indi- 
cates the Mott-like regime, since otherwise it would not 
be possible to switch back (almost) adiabatically to the 
initial state. Thus, this figure provides a glimpse at a 
quantum phase transition induced by "dressing" a mat- 
ter wave, although, of course, a truly sharp "transition" 
cannot be achieved with N = M = 7. 

But the small system already is sufficently rich to 
demonstrate that an ideal outcome cannot be taken for 
granted: Fig.[2]shows a momentum distribution obtained 
in the same manner for a lower frequency, Tilo/J — 12. 
Whereas a signature like this might not be distinguish- 
able from that in Fig. [1] in an actual experiment for times 
up to i2, here the initial pattern is not restored at ts, 
indicating severe deviations from the desired adiabatic 
following. 

Inspection of the corresponding quasienergy spectra 
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FIG. 3: Part of the quasienergy spectrum for the almost ideal 
case considered in Fig. [1] The arrow marks the quasienergy 
of the Floquet state evolving adiabatically from the ground 
state of the undriven Bose-Hubbard system. Computed with 
N = M = 5. 




FIG. 4: Part of the quasienergy spectrum (A'^ = M = 5) for 
the thwarted case considered in Fig. [21 the arrow indicates 
the quasienergy associated with the ground state. Observe 
that here there are several active resonances, corresponding 
to pronounced avoided level crossings. 



immediately reveals the reason for the different dynam- 
ics found in both cases. In Fig. [3] we depict a part of 
the quasienergy spectrum for the frequency employed in 
Fig. [U whereas Fig. |4] shows the spectrum for the lower 
frequency underlying Fig. [21 both spectra have been com- 
puted with iV = M = 5 by solving the eigenvalue equa- 
tion l[?71) . In the first case, the quasienergy line emanat- 
ing from the ground state is not visibly affected by other 
states, indicating the absence of active resonances and 
thus enabeling the adiabatic return observed in Fig. [2 
essentially; only one single instantaneous Floquet state is 



populated during the entire process. In contrast, in the 
second case the quasienergy level originating from the 
ground state undergoes several large avoided crossings. 
Incomplete Landau-Zener transitions at these avoided 
quasienergy crossings then lead to a significant popula- 
tion of the anticrossing Floquet states [22i | , rendering an 
adiabatic return to the initial state impossible. 



VI. CONCLUSIONS 

Ultracold atoms in time-periodically modulated opti- 
cal lattices give rise to "dressed matter waves" , in anal- 
ogy to the dressed atoms known from atomic physics [l| . 
Such dressed systems acquire properties quite different 
from their "bare" antecedents, the modification of atomic 
Lande g-factors setting a prominent example An ef- 
fect closely related to this g-factor modification is a tran- 
sition of ultracold bosonic atoms in a modulated optical 
lattice from a superfluid to a Mott-insulator-like state in 
response to a variation of the modulation strength; this 
transition is mediated by a modification of the nearest- 
neighbour hopping matrix elements relying on precisely 
the same mechanism as that of the g-factors. An experi- 
mental verification of this proposal involves adiabatic fol- 
lowing of the many-body Floquet state originating from 
the ground state of the bare system; such adiabatic fol- 
lowing is a quite tricky concept in the context of driven 
matter waves. The question to what extent multiphoton- 
like resonances can be avoided (or perhaps deliberately 
be induced and exploited) is open to experimental inves- 
tigation. 

The long-term perspective of these considerations, 
however, seems to lie elsewhere. Just as the g-factor 
modification is but one facet of the dressed-atom picture, 
there are further possibilites of controlling the state of a 
matter wave in an optical lattice by time-periodic forc- 
ing. Our present scheme defines a first cornerstone; if 
achieved, more demanding ones can follow. In particu- 
lar, it might be interesting to resonantly couple different 
Wannier states located at the same site, and thus to open 
up new ways of quantum state engineering. The experi- 
mentally established fact that Bose-Einstein condensates 
in optical lattices can be subjected to strong forcing with- 
out destroying their phase coherence [J] is a sound cause 
for optimism. 
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